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Abstract 

This study proposes an algorithm to synthesize controllers for the power management on board 
hybrid vehicles that allows the vehicle to reach its maximum range along a given route. The algorithm 
stems from a level-set approach that computes the reachable set of the system, i.e., the collection 
of states reachable from a certain initial condition via the computation of the value function of 
an optimal control problem. The discrete-time vehicle model is one of a particular class of hybrid 
vehicles, namely, range extender electric vehicles (REEV). This kind of hybridization departures 
from a full electric vehicle that has an additional module - the range extender (RE) - as an extra 
energy source in addition to its main energy source - a high voltage battery. As an important 
feature, our model allows for the switching on and off of the range extender and includes a decision 
lag constraint, i.e., imposes two consecutive switches to be separated by a positive time interval. The 
approach consists in the introduction of an adequate optimal control problem with lag constraints 
on the switch control whose value function allows a characterization of the reachable set. The 
value function is in turn characterized by a dynamic programming algorithm. This algorithm is 
implemented and some numerical examples are presented. 

Keywords: Optimal control, Range Extender, Hybrid vehicles, Reachability analysis. 



1 Introduction 

Electrified automobile power chain technology faces a rapid development. The interest for this technology 
arises as an increasing number of automakers wish to adopt vehicles with an electrification of the power 
chain as a viable solution for reducing greenhouse gas emissions worldwide to meet stringer regulative 
legislation and consumers' demand. There is an effort of major constructors to deploy fully electric 
vehicles (EVs) as early as possible in some car market segments. This work focuses on a specific class 
of hybrid electric vehicles, namely, range extender electric vehicles (REEV). The range extender electric 
vehicle model consists of a small dimensioned internal combustion engine (ICE) as an extra energy source 
acting as a range extender in addition to the main energy source, namely, a DC battery. The two energy 
sources are assembled as in a series hybrid electric vehicle, which means that the ICE is not connected 
mechanically to the drive shaft. Instead, a generator transforms the mechanical energy produced in the 
ICE into electric energy that can be directed towards the electric motor or the battery. Because of the 
relative low power of the ICE, this configuration does not allow the vehicle to rely solely on the range 
extender and thus the vehicle cannot assure traction if there is no energy in the battery. Although our 
controllers have to take into account the particularities of this class of hybrid vehicles, the formulation 
developed in this study is general enough and may be adapted to other architectures, including series 
and parallel hybrid electric vehicles. 

The goal of this work is compute a way of operating the range extender to enable the vehicle to 
reach its maximum range. The REEV is modeled as a discrete-time dynamical system in which the 
state vector represents the energy capacities of the two different energy sources. In this setting, the 
study aims at finding the control sequence of the two energy sources that allows the vehicle to reach the 
furthest possible point of a given route or to reach the driver selected energy state. 

To the extent of our knowledge, there are no studies considering the driving range of the vehicle 
as an optimization criterion. Usually, in the context of power management strategies the optimality 
of a strategy is defined in terms of the total fuel consumption during a driving cycle. We refer here 
to [TJ|2j[5j[6j|9] . Different techniques are used in each of these works, as well as different assumption on 
whether one has access to a priori driving information. A common feature of the works above is the 
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fact that the power management synthesis uses a predetermined driving range, typically given as one 
considers a particular driving cycle. 

The main idea in our work is to have no a priori information about the vehicle driving range. A 
typical setting is when the driving cycle is too long to be completed with the current energy on-board 
the vehicle. In the case where no power management strategy can ensure the cycle completion, the 
utilization of a fuel consumption criteria does not seem to lead to well-posed posed optimal control 
problem. For instance, dynamic programming-based algorithms cannot be implemented, since they rely 
on a backward computation of the value function starting from the "last" reachable point on the route 
which may be unknown. This work proposes a method to synthesize a power management strategy that 
maximizes the driving range on a cycle supposedly too long to be completed. 

This paper is organized as follows: firstly, it describes the REEV model motivating the study and 
states the associated hybrid optimal control problem. Then, the reachable set and the value function 
are defined and a dynamic programming principle for the value function is obtained. It follows that the 
dynamic programming principle is used to derive an dynamic programming algorithm to compute the 
value function thus characterizing the reachable set of the system. Lastly results of numerical simulations 
evaluating the autonomy both of a REEV toy and a realistic model are presented. 

2 Problem Settings 
2.1 Vehicle Model 

This section discusses the model of the power management system of a REEV that is to be optimized and 
introduces the notation used throughout the document. A REEV is a vehicle that combines a primary 
power source - a HV battery - and a small dimensioned (powerwise) secondary power source - in our case 
an ICE. The traction (or propulsion) of the vehicle is performed by an electric motor connected to the 
vehicle's wheels through a reduction gear. Both energy sources can supply some share of the necessary 
power to meet the driver demand. Additionally since the model considers a range-extender electric 
vehicle type, it cannot rely solely on the RE power to drive the vehicle. Thus, the model considers that 
the vehicle's traction capability is conditioned to the existence of some electric energy in the battery. 

The architecture is that of a series hybrid electric vehicle, which means that the ICE is not mechan- 
ically connected to the transmission. Instead, a generator transforms the mechanical energy produced 
in the ICE into electric current that can be directed towards the electric motor or charge the battery. 

The ICE state is controlled by a discrete sequence of switching orders. An important feature of the 
REEV model is a time lag 5 > imposed between two consecutive decisions times. From the physical 
viewpoint, this assumption incorporates the fact that frequent switching of the RE is undesirable in order 
to avoid mechanical wear off and acoustic nuisance for the driver. This lag has an order of magnitude 
of about 120s. 

The controlled variables are the state of charge (SOC) of the battery, the remaining fuel on the ICE 
reservoir and the ICE state on or off. 

This setting considers that the optimal power management strategy is synthesized for a specific trip. 
In this direction, we assume that information about the vehicle's future power demand is known. In 
particular, it is assumed that the vehicle future speed, i.e. the speed profile of the vehicle, as well 
as the relief of the route followed by the vehicle are known. As a consequence, the problem can be 
stated in a deterministic setting. This may seem a very strong assumption. However, range extender 
electric vehicles are expected to be equipped with embedded navigation systems (NAV) that actually 
have information about expected speeds and topological data. This information can be exploited in 
order to construct an expected speed profile with a certain degree of confidence. 

In the following of the paper, index k = 1, • • • , k shall index a position in the route. Since the vehicle 
speed is assumed to be known, on can easily obtain, from a known location k, the time to attain it. 
More precisely, denote by dk the distance between points k—1 and k and o~k > the (constant) vehicle 
speed between points k—1 and k. The total time to drive between k — 1 and k is then simply dfc/ofc. 

The state vector is composed of the vehicle energy state and the ICE state. The vehicle energy 
state is a two-dimensional vector y e X = IR 2 , where y — (2/1,2/2) denotes the state of charge of 
the battery and the fuel available in the range extender module respectively. Each of these quantities 
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are the image of the remaining energy in the battery and in the RE respectively. Moreover, because 
energy capacities of the battery and the fuel tank are limited (and normalized between 1 - full reservoir 
- and - empty reservoir), the set of admissible states is defined as K = [0, l] 2 . Remark that we 
define the state space as all the plane X = IR 2 and introduce a set of admissible states K. The set 
K represents the "physically admissible" SOC and fuel levels, considering that they are normalized 
capacities. This is a modeling choice that later enables us to deal easily with the state constraints in the 
dynamic programming principle (in opposition to imposing "hard" constraints on the state space like 
X = [0, l] 2 ). The RE state is denoted by q G Q = {0, 1} and indicates whether the RE is off (q = 0) or 
on (q = 1). 

The control variables are the output power of the ICE and the order to switch the ICE's state (on or 
off). The first is referred to as the continuous control and the latter as the discrete (or switch) control. 
The continuous control is applied at all points k = 0, • • • , k — 1 of the route and is denoted Uk G U{qk). 
The idea is that at each point k, also called the decision stages, the controller is able to control the 
power split of the vehicle. Indeed, in a series architecture the power split is controlled only by setting 
the output power of the range extender. That is because, for a given required power at the wheel P*, 
imposed by the driver, once the range extender output power u is chosen the remaining required power 
P* — u is automatically supplied by the battery. 

The discrete control w is a sequence of switching decisions w = {(toi, Si), • • • , (iojv, sat)}, where for 
all j — 1, • • • ,N, Si G {1, • • • , k — 1} and wj G W(q Sj ) C Q. The sequence of switch controls {wj}jLi 
(designating the new mode of operation of the ICE) is associated with the sequence of switching positions 
{s^jLi so that each decision Wj is exerted at node Sj. The set W(q) is a set value function of the actual 
discrete state and contains all possible new discrete states. In the vehicle application, where there are 
only two discrete modes of operation (on or off), this set is simply W(l) = {0} and W(0) = {1}. The 
switching decision is introduced so one is able to stop the range extender whenever there are large idle 
periods, thus saving fuel. More precisely, if the ICE is running and no output power is required, there 
is still an idle fuel consumption, as energy is required to keep the engine running to avoid stalling. In 
the case of extended idle periods, it may be worth to stop the ICE and pay a fuel overhead to restart 
the engine whenever needed. 

The switch lag condition imposes that two consecutive switch decisions must be separated by a time 
interval of 8 > 0, i.e. 

s j+i , 

^>5. 



E 



1 ai 



Given a discrete control sequence to and a discrete state qj, the continuous control steers the state 
according to 

Uk-i = fk(Vk,Uk-i,qj) 
and at isolated points {sj}, the discrete state changes according to 

'/.. i '/•:</...• 

These expressions mean that: given a discrete control sequence w, let qj be a mode of operation activated 
at some point Sj + \. Then, the dynamic j '.(•,-, qj) operates on the state between points k = Sj+i, ■ ■ ■ , Sj 
and the trajectory of the discrete state is simply = qj for k between Sj+i and Sj. When point Sj is 
reached, the discrete dynamics undergoes a switch associated with the switch control Wj and changes to 
qj-i — g(qj,Wj), activating the new mode of operation given by /.(•, -,qj-\) . Observe moreover that 
the process is modeled and described in backward fashion. 

At this point, wc introduce the notion of hybrid control and admissible hybrid control. The hy- 
brid control is introduced for notational convenience. It is a control regrouping both the range ex- 
tender power u and the switch decision w. More in detail, given a discrete control sequence w = 
{(wi, Si), • • • , (u>n, sn)} and a continuous control sequence u = (u\, • • • , u K -i), the hybrid control is de- 
noted by a — (u,w) and groups the continuous and discrete controls. Denote by A the set of all hybrid 
controls. As mentioned, the hybrid control a must contain controls (u, w) such that Uk G U(qk) and 
Wj G W(q Sj ) for j = 1, • • • , N. However, the admissibility of each discrete and continuous control does 
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not imply the admissibility of the hybrid control without additional structural conditions. By structural 
conditions we refer to the fact that the decision nodes Sj must be increasing with j and that the asso- 
ciated decision times must respect the lag condition. The structure imposed on the admissible hybrid 
controls a <E yl adm allows the system to engender admissible trajectories. An admissible trajectory can 
be regarded as a "physically" admissible trajectory. Such trajectories correspond to the evolution of the 
SOC and fuel levels whenever the system is steered by an admissible control sequence. However, the 
hybrid control admissibility condition is not well adapted to a dynamic programming principle formu- 
lation, needed later on. In order to include the admissibility condition in the optimal control problem 
in a more suitable form, we introduce a new state variable tt. Recall that the decision lag condition 
implies that new switch orders are not available until a time S after the latest switch. The new variable 
is constructed such that, at a given position k, the value tt^ measures the time since the last switch. 
The idea is to translate the structural conditions into conditions satisfied by this new variable and treat 
them as state constraints. Thus, if TTk < $ all switch decisions are blocked and if, conversely, > S the 
system is free to switch. For that reason, this variable can be seen as a switch lock. Now, given n > 0, 
k = 1, • • • , K, a discrete control sequence w, final conditions x G X and q € Q, the hybrid dynamical 
system is written as 

Vk-i = fk(yk,u k -i,qj), k = k, Vk = x (1) 

Qj-i = g(qj,ivj), i = i)--- >-W, qN = q- (2) 

( <* + £•=! if k< Sl 

= ^ = inf * - ^ if k > Sl (3) 

Indeed, once the discrete control is given, the trajectory irj. can be determined. Proceeding with the 
idea of adapting the admissibility condition we wish to consider ir K — p, with p G P := (0, T], where 
T is the total travel time, the final value of the switch lock variable trajectory and impose the delay 
condition under the form 7r s ._x > S for all Sj (notice that n s . = by construction). Proceeding as 
such, these conditions suffice to define an admissible discrete control set. So, while optimizing with 
respect to admissible functions, one needs only look within the set of hybrid controls that engender a 
trajectory ir^ with an appropriate structure. Denote the solutions of ([l])-([3]) with final conditions x,q,p 
by Vx,q,p;K,qx,q,p;K and TT x ,q,p;K- Given k>0, ?eQ and peP, define the admissible trajectory 

Set ^(0,k) aS 

S (o!k) = {y*,q,p;K I a = ( u A w j> s j}f=i) S A, y x ,q,p-K solution of @-([2]), 

^x,q,p; K Solution Of 7T K = p, 7T Sj _l > 6, j = 1, • • • , N}. (4) 

2.2 Reachability Problem and Autonomy 

Let Xq C X be the set of allowed initial states, i.e. the set of states from which the system ([lj-([2| is 
allowed to start. Define the reachable set as the set of all points attainable by y after driving until point 
k starting within the set of allowed initial states Xq to be 

Rf = {x I G Q x P, y x ^ p . k G S***, 

(yx,q, P ;k)o G X°, and (y x , qiP] k)e G K, 9 = 0, • • • , k] (5) 

In other words, the reachable set R^° contains the values of the final energy state (y x ,q;k)kt regardless 
of q,p at point k, for all admissible trajectories - i.e., trajectories obtained through an admissible hybrid 
control - starting within the set of possible initial states X that never leave set K. In particular, the 
information contained in ^ allows one to determine the first point where the reachable set is empty. 
More precisely, given X Q C X, define k* > to be 

k* = inf{fe > | Rf c 0}. (6) 

The position §6§ is identified as the autonomy of the hybrid system ([l])-([2]). Indeed, one can readily see 
that if no more admissible energy states are attainable after k* , any admissible trajectory must leave 



REACHABILITY OF DELAYED HYBRID SYSTEMS USING LEVEL-SET METHODS 



5 



set K of admissible values beyond this time. Therefore, k* is interpreted as the last point in the path 
in which there is some usable energy on-board the vehicle. 

In the case of range extender electric vehicles, the vehicle stops operating in nominal mode when the 
battery is depleted, even if there is some fuel left in the tank. If the driving conditions allows all the 
fuel to be consumed, then the last reachable energy state of the vehicle is the state (0, 0), where there is 
no charge in the battery and no fuel in the tank. However, this is by all means the sole last reachable 
energetic state in all situations. Indeed, because the range extender is power-wise small dimensioned, 
it may be the case that the engine cannot consume energy quickly enough to empty the fuel tank. For 
instance, if the vehicle is driving through a steep slope, the battery may be depleted before all the fuel 
is consumed resulting thus in a final state (0,2/2), 2/2 > 0. The reachable set R£ is associated with a 
minimum time function T : X — > ]R + U {00}, defined as 

T(x) = inf{fe > I x e Rf } (7) 

The minimum time function associated to a point x in the state space is the minimum time that any 
admissible trajectory needs to reach x departing from set Xq while respecting the state constraints K. 



2.3 Optimal Control Problem and Dynamic Programming Principle 

In order to characterize the reachable set R*° this paper follows the classic level-set approach 7|. The 
idea is to describe ^ as the negative region of a function v. It is well known that the function v can be 
defined as the value function of some optimal control problem. In the case of system Q-Q, v happens 
to be the value function of a (discrete time) hybrid optimal control problem. Consider a Lipschitz 
continuous function (f> : X — > IR, such that 

<f>(x) < o x e X°. (8) 

Such a function always exists - for instance, the signed distance function dx from the set X°. 

For a given point k > and hybrid state vector (x, q,p) G X x Q x P, define the value function to be 

v° k (x,q 7 p)= mi p {<f>((y x , q , p .,k)o) I {Ux,q, P -k)e € K, V0 £ {0, • ■ • ,k}}. (9) 

Observe that ^ works as a level-set to the negative part of Indeed, since ^ contains only admissible 
trajectories that remain in K, by (JsJ) implies that v®(x,q,p) is negative if and only if (y x ,q,p-.k)o is inside 
Xq, which in turn implies that x E R k . 

Now we proceed to describe the state constraints in the similar manner. The set X\K of state 
constraints is treated as an obstacle in the state space, i.e., the state trajectory cannot pass through the 
obstacle set K. We recall that K = [0, l] 2 is the normalized SOC and fuel quantities. More precisely, 
define a Lipschitz continuous function ip : X — > TR to be 

if{x) <0^xeK. (10) 

Then, for a given k > and (x, q,p) G X x Q x P, define a total penalization function to be 

Jk(x,q,p;y) = ( <t>((y x ,q, P ;k)o) \J max v>{{y x , q ,p-k)e) 

\ #£{(),••' ,k\ j 

and then, the optimal value : 

Vk{x,q,p)= inf J k (x,q,p;y). (11) 

Here and in the rest of the paper, we denote a V b := max(a, b). 

The next proposition certifies that ^ is indeed a level-set of ([9| and (11). 



Proposition 2.1. Let Xq C X. Define functions 4> et ip respectively by ( |10| and fl8j). Define value 
functions v° et v respectively by (11) and Then, for k > 0, the reachable set ^ is given by 

R*° = {x I 3{q,p) e Q x P, v k (x,q,p) < 0} - {x \ 3(q,p) e Q x P v° k (x,q,p) < 0} . (12) 
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Proposition (2.1 1 sets the equivalence between ^ and the negative regions ol Q and (11). In 
particular, it states that it suffices to computes v or vq in order to obtain information about R Xa . In 
this sense, this paper focuses on ( 11 ), which is associated with an optimal control problem with no state 
constraints. The value function can also be used to characterize the minimum time function. In addition 
to the minimum function, the next proposition introduces the extended minimum time function, which 
takes into account all state variables x,q,p. This extended function is needed later to the controller 
synthesis. 



Proposition 2.2. The minimum time function T : X — > IR + U {oo} is given by 

T(x) = inf{k | 3(q,p) G Q x P, v k (x,q,p) < 0} . 
The extended minimum time function T 1 : X x Q x P — > ]R + U {oo} is given by 

T'{x,q,p) = inf {k | v k (x,q,p) < 0} . 



(13) 
(14) 



Observe that equation (13) is an immediate consequence of the level set identities in proposition 2.1 



Relation ( 14 ) characterizing the extended minimum time function using the value function is derived in 
a similar fashion. 



Before stating a dynamic programming principle satisfied by (11), we introduce some preliminary 
notations. Given k > 0, set JC = {1, • • • , n}, H = XxQxPxK, and denote its closure by fi. (We 
recall that, concerning the vehicle application, X = IR 2 , Q = {0, 1},P = (0, T], where T is taken to be 
the total travel time). As such, |TT] ) is defined as 

v : H -> IR 

(x,q,p,k) v k (x,q,p) 

Define 

V(fi) := {v | v : O -> IR, v bounded }. 
Now, define the non-local switch operator M : V(f2) — > V(f2) to be (adapted from the ideas in [10] ) 

(Mv k )(x,q,p) = inf v k (x,g(w,q),p') 

wGW(q) 
p'>5 

The action of this operator on the value function represents a switch that respects the delay constraint. 
Thus, they only make sense when p — 0. The next proposition states the dynamic programming principle 
verified by (ITT). 



Proposition 2.3. The value function (11) satisfies the following dynamic programming principle: 
(i) For k=0, 



(ii) For p = , 
(in) Otherwise, 



v (x,q,p) = 4>{x) \/<p(x), V(x,q,p) G X x Q x P, 
v k {x,q,0) = (Mv k )(x,q,0), V(x,q,k) G X X Q X AC, 



Vk(x,q,p)= inf <v k - 1 [ f k (x.u k - 1 ,q),q,p ■ 



dk 



(15) 
(16) 

(17) 



The dynamic programming principle allows one to compute the value function and consequently 
characterize the reachable set R k ° ■ Once the value function is obtained, one can synthesize a controller 
using the minimum time function T . The computation of the value function and the subsequent controller 
synthesis are obtained in a classic fashion by a dynamic program ming algorithm. The algorithm is the 
projection of the dynamic programming principle (cf. proposition 2.3 ) over a grid X* x Q x P# of the 
discretized state space. The dynamic programming algorithm outputs a controller (u* , w*) that controls 
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the range extender power output and its mode of operation to reach the destination point after a time k 
such that the final state is (x',q',p') departing from a state in Xq. Additionally, the trajectory remains 
inside set K, which in our case amounts to respect the fuel and SOC capacities. One can observe that 
the extended minimum time function plays a major role in the synthesis of the controller. Indeed, the 
value function numerical values have no physical meaning since the obstacle and target functions <f> and 
ip are arbitrary. Only the sign of the value function is needed to determine the reachable set. Instead, 
the physical information is encoded in the minimum time function. 



3 Numerical Simulations 



Proposition 2.3 provides the expressions that allow the computation through a dynamic programming 



algorithm. Once the value function obtained, one can use ( 12 ) to characterize the reachable set and 
^ to obtain the autonomy of the system. First, a simple vehicle model is used in the simulations. In 
addition, a constant vehicle speed is assumed during this first simulation phase. This model allows an 
analytic evaluation of the autonomy and is suitable for an a posteriori verification of the results. After 
the algorithm is tested using this simple configuration, we introduce a realistic vehicle model and a 
realistic speed profile. 

In this simple setting, we assume that the vehicle follows a constant speed of a = lm/s and the 
distance between any two points k—1 and k is dk — lm. We consider that X = 1R 2 and Q = {0, 1}. The 
switch dynamics are given simply by g(w,q) = w ensuring that w imposes the new mode of operation 
at every switching time. The continuous state dynamic models a simple REEV charge behavior. The 
evolution of the energetic state (x, y) £ IR 2 - respectively, the battery SOC and the range extender fuel 
- is given by f(u,q) = (— a x + qu,—q(a y + u)), where a x ,a y > are constant depletion rates of the 
battery electric energy and the reservoir fuel (whenever the RE is on), respectively. The coefficient a x 
models a non-negative power requirement throughout the driving cycle and therefore, a non-negative 
SOC consumption. However, the net SOC variation can be made zero or even positive (i.e., charge 
the battery) by the action of the output power of the range extender u. To do so, the range extender 
must be running and fuel must be consumed. The factor a y represents the fuel consumed to generate the 
necessary power to counter the inertia of moving parts of the range extender while u is the effective share 
of power transfered to the electric motor. In order to simplify the analysis we consider the conversion 
factors from units of power to units of (normalized) fuel mass and percentage of battery SOC to be 
unitary. 

The engine maximum power is denoted by u max and the minimum power is set to be 0. 

Considering this simple dynamics, an exact autonomy of the system can be evaluated analytically. 
Since there are no non-linearities in the model, an optimal strategy to reach the furthest point is to empty 
completely the fuel as soon as possible, thus obtaining the maximum amount of energy available. Given 
initial conditions (xo, yo) the shortest time to empty the fuel reservoir is given by t* = yo/(a y + M max ). 
The SOC evaluated at this instant is given by x(t*) = x(0) — t*(a x — u max ). If x{t*) < 0, it means the 
fuel cannot be consumed fast enough before the battery is depleted. This condition can be expressed 
in terms of the parameters of the model as xo(a y + w max ) < yo(&x — u max ). In this case, the autonomy 
is given by T° = — ^ — . Time T° is the time needed to spend all the energy on board, using the 
range extender at full power. In the other hand, if x(t*) > 0, all the fuel in the range extender can be 
consumed and there is still some time until the battery is emptied. In this case, the autonomy of the 
system given by T 1 — a:o+ ^ maxt . 

Figure [3] illustrates a controller synthesized for a given initial condition. One can readily see that the 
reachable set, and in particular, the trajectory remain inside the authorized region of the state space 
K = [0, l] 2 . Also, the controllers respect the delay condition 8 = 1 between two consecutive switches. 
The initial set is taken to be a ball around (xo,yo) m each case for numerical reasons. 

Table [3] presents the error between the theoretical and computed autonomy as well as the CPU time 
required for two initial conditions in the instance using a x = 0.10, a y = 0.15, w max = 0.07 and 6 = 1. 
One can observe the convergence of the algorithm towards the theoretical value as Ax approaches 0. 
Empirically, a negligible sensibility with respect to Ap is observed. 

'-Intel Xeon E5504 @ 2 x 2.00GHz, 2.99Gb RAM. 
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SOC time (s) 

Figure 1: Figure showing the trajectory of the energy state from a valid initial state until the last 
reachable point. The plot on the left groups the contour levels of the minimum time function, the 
obstacle and initial sets. On the right, the figure displays the range extender power control u and 
the switch sequence w controlling the ICE state. The synthesis is obtained using step sizes of Aa: = 
0.0245, Ap = 0.4211, At = 0.4. 





Ax 


£ 


CPU 1 running time(s) 




0.05 


3.291 


2.48 




0.04 


1.154 


2.89 


(0.5,0.5) 


0.03 


0.218 


4.21 




0.02 


0.081 


8.55 



Table 1: Convergence results and running times. Fixed At — 0.4 and Ap — 0.5. Instance considering 
a x = 0.10, a y = 0.15, u max = 0.07 
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Initial SOC, fuel 


(0.3,0.3) 


(0.6,0.4) 


Max range 
EV range 

Relative range increase 


45.126 km 
22.045 km 
104.70% 


78.405 km 
48.209 km 
62.64% 


Range extender op. cost 


11.70 €/100km 


11.92 €/100km 



Table 2: Results for the autonomy of the range extender electric vehicle. The table presents the maximum 
range of the vehicle when equipped with a range extender (Max range) and the range of the same vehicle 
when operating in a purely electric mode (EV range). The financial cost of operation of the range 
extender is estimated, based on a price of 1.5 €/l. 

These remarks about the behavior of the algorithm allow us to proceed to the study of the reachable 
set in the realistic vehicle model. The procedure for implementing the simulations is as follows. We 
consider a driving profile consisting of n > nodes and k — 1 links. At each node fc, entry to the link k 
there is an associated link speed and the length of the link d^. Once cr^jdj^ are known, one can fix 
the time steps Atk used in the algorithm. Notice that the time steps depend on the link, the problem 
being non-autonomous in this case. 

We proceed to the computation of the autonomy using a realistic vehicle model. Table [3] summarizes 
the main results. Figure [3] regroups the minimum time function, the synthesized controller that allows 
the vehicle to reach its maximum range and the corresponding evolution of the SOC and fuel. The 
autonomies computed in both cases can be compared to the autonomy of the vehicle under a purely 
electric utilization. In the first case, the utilization of the range extender roughly doubles the range of 
the vehicle (a 104.7% range increase). When considering the initial condition with a SOC at 60% and 
fuel at 40%, the range extender increases the vehicle range by over 60%. We also consider the total 
operating financial cost (given in euros in this study) of the range extender. This cost is based on an 
estimated r\ = 1.5 €/l of fuel. The range extender operating cost is then normalized with respect to 
100km to yield a more meaningful result. We stress that we take into consideration in the cost of the 
range extender operation only the additional distance due the range extender. Therefore, the values 
of operating cost should be read as euros per 100km of range extender utilization instead of 100km of 
driving distance. In our model the reservoir has a total capacity of 61. 

4 Conclusion 

In this study we present an algorithm to determine the autonomy of a hybrid vehicle. The algorithm 
uses information about the driving profile of the vehicle, obtained typically through a navigation sys- 
tem. Once the destination is fixed and if there is not enough energy on-board the vehicle to reach the 
destination point, the algorithms computes the longest driving distance and synthesize the associated 
control sequence of the controllable power source. In our application, we focus on the range extender 
electric vehicle-class of hybrid vehicles. Our model considers an internal combustion engine as a range 
extender and a high voltage DC battery as main energy source. The dynamic evolution of the state is 
considered to be a discrete-time evolution. 

The algorithm computes the reachable set, set of reachable states, given a final time and initial 
state, via a level-set approach. The reachable set is characterized as the negative region of a value 
function associated to an optimal control problem. To compute the value function, we formulate an 
optimal control problem with obstacle and we state a dynamic programming principle verified by the 
optimal control problem. Then, the value function can be computed using a (deterministic) dynamic 
programming principle. 

Numerical simulations of this procedure are performed first using a hybrid vehicle toy model and 
then using a realistic range extender electric vehicle model. The results shows that, depending on the 
initial conditions, operating the range extender as described in this paper may yield a 100% increase in 
the vehicle range, with respect to a purely electric vehicle. 
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Figure 2: Figure showing the trajectory of the energy state from an valid initial state until the last 
reachable point. The plot on the left groups the contour levels of the minimum time function, the 
obstacle and initial sets. On the right, the figure displays the range extender power control u and the 
switch sequence w controlling the ICE state. The synthesis is obtained using step sizes of Ax = 0.0308. 
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